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Abstract. The explicit semiclassical treatment of logarithmic perturbation theory 
for the bound-state problem within the framework of the radial Klein-Gordon equation 
with attractive real-analytic screened Coulomb potentials, contained time-component 
of a Lorentz four-vector and a Lorentz-scalar term, is developed. Based upon 
?i-expansions and suitable quantization conditions a new procedure for deriving 
perturbation expansions is offered. Avoiding disadvantages of the standard approach, 
new handy recursion formulae with the same simple form both for ground and excited 
states have been obtained. As an example, the perturbation expansions for the energy 
eigenvalues for the Hulthen potential containing the vector part as well as the scalar 
component are considered. 
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1. Introduction 

Static screened Coulomb potentials have been widely used in nuclear and particle 
physics, atomic physics, solid-state physics and chemical physics. Usually the bound- 
state problem with these potentials is considered within a non-relativistic framework. 
Nevertheless, relativistic effects for a particle under the action of such potential could 
become important, especially for strong coupling. 

Several attempts have been made to describe relativistic systems in a central field 
due to a time-component of Lorentz four-vector and a Lorentz-scalar interaction within 
the framework of the Klein-Gordon equation. However, for almost all potentials this 
equation is not exactly solvable which compels to resort to some approximation methods. 

A number of such approaches to solving the radial Klein-Gordon equation in 
analytical expressions have been developed, including, in particular, the use of the 
relativistic hypervirial and Hellmann-Feynman theorems the 1/A^-expansions [2- 
9], the method of Regge trajectories [10-12], the elements of an 5*0(2,1) Lie algebra 
||13| , |I4[ and a perturbation scheme based on a comparison equation [|1^ . 

Despite such a variety of methods one of the most popular techniques is still 
logarithmic perturbation theory [16-18]. Within the framework of this theory, the 
conventional way to solve the quantum-mechanical bound-state problem consists in 
changing from the wave function to its logarithmic derivative and converting the Klein- 
Gordon equation into the nonlinear Riccati equation. In the case of ground states, the 
consequent expansion in a small parameter leads to a hierarchy of simple equations 
that permit us to derive easily the corrections to the energy and wave function for each 
order. However, when radially excited states are considered, the standard approach 
Tsj becomes extremely cumbersome and, practically, inapplicable for describing higher 



orders of expansions. At the same time the evaluation of perturbative terms of large 
orders is needed for applying modern summation procedures because the obtained series 
are typically divergent. 

The above mentioned weakness of the standard approach is caused by factoring 
out zeros of the wave functions with taking into account corrections to the positions 
of the nodes. On the other hand, the radial quantum number, n, that is equal to the 
number of nodes of the wave function, most conveniently and naturally is introduced in 
the consideration by means of quantization conditions as well as in the WKB-approach 



||T9| , pO| . However, since the WKB-approximation is more suitable for obtaining energy 
eigenvalues in the limiting case of large quantum numbers but perturbation theory, on 
the contrary, deals with low-lying levels, the WKB quantization conditions need change. 

Recently, a new technique based on a specific quantization conditions has been 
proposed to get the perturbation series via the semiclassical ^-expansions within the 
one dimensional Schrodinger equation |^T|, 

In this paper, we would like to extend similar technique to the bound-state 
problem for the Klein-Gordon equation with the Coulomb-type potential, contained 
time component of a Lorentz four- vector and/or a Lorentz scalar term, that receives 
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attention as possible model of quark confinement. 

There are several advantages to described approach: no wave functions or matrix 
elements need to be previously calculated and derived recursion formulae for obtaining 
perturbation corrections have the same simple form both for ground and excited states. 

The proposed technique can be regarded as a further investigation of a part assigned 
to a rule of achieving a classical limit for radial and orbital quantum numbers in the 
construction of the semiclassical methods for the Klein-Gordon equation. If the WKB- 
approach is realized using the rule, h ^ 0, n oo, I ^ oo, hn = const, M = const, 
application of the condition h ^ 0, n = const, I — > oo, hn 0, hi = const was proved 
to lead to the method of 1/A^-expansion p|, |^ and the method of Regge trajectories 



111 , [121 . Here we address ourselves to the alternative possibility: h 0, n = const, 
I = const, hn 0, hi ^ 0, that results in the explicit semiclassical treatment of the 
logarithmic perturbation theory for the radial Klein-Gordon equation. 

This paper is organized as follows: section II contains a general discussion of 
the method and necessary assumptions in semiclassical treatment of the logarithmic 
perturbation theory. In section III quantization conditions obtained are used for deriving 
the recursion relations for perturbation expansions. Section IV demonstrates that 
described approach restores the exact results in the case of the Coulomb potential 
and gives the example of its explicit application to the bound-state problem in field 
of the Hulthen potential. This potential, apart from its initial interest as a possible 
form of nuclear interaction, is exactly solvable for the / = states, thereby providing a 
consistency check for our perturbation technique. Finally we present a brief summary. 

2. The method 

In this section, we study the bound state problem for a relativistic scalar particle 
moving in the field of an attractive central real-analytic screened Coulomb potential, 
admitted bounded eigenf unctions, having in consequence a discrete energy spectrum. 
This potential has a Coulomb- like behaviour at the origin, caused by the time component 
of a Lorentz four-vector, V{r), and/or a Lorentz-scalar term, W{r), which in general 
can be written as 

1 oo 1 oo 

V{r) = -Y^ V, r' Wir) = W, / . (1) 

Notice, that in the last case, when the Coulombic singularity is contained only in one 
part of potential, another part may be a smooth function described by equation with 
Vo = OotWo = 0. 

In what follows, a scalar potential will be included in the mass term m{r), by 



analogy with "dynamical mass" models of quark confinement p3|, P4 



W(r) 

m{r) =m + . (2) 
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Then the reduced radial part of the Klein-Gordon equation takes the form 



1 



h^R"ir) = ^m\r)c' - ^[E - ^(r)]^ + + | R{r) (3) 

where the particle wave function is \E'(r) = R{r)/r. 

As is customary in the logarithmic perturbation theory, we apply the substitution, 
C{r) = hR'{r)/R{r), and go over from (^) to the Riccati equation 

hC'ir) + C\r) = !^!iii±i) + rn\ry - \\E - V{r)]^ . (4) 

Within the framework of the standard approach to the logarithmic perturbation 
theory, all quantities of the Riccati equation are expanded in formal series in powers of 
a small screening parameter A, incoming in a potential functions of the form V{r) = 
r~^/(Ar). After such an expansion the screening parameter appears in consideration 
in combination with derivatives of screening function as A* d*//dr*. Taking into account 
this fact, in our approach we do not single out the screening parameter, but incorporate 
it into coupling constants Vi and Wi from (||). 

Besides, because the Riccati equation (^) has in the relativistic case the same 
structure as in the nonrelativistic one |21|, these coupling constants, Vi and Wi, 
appear in common with powers of Plank's constant. Therefore the perturbation series 
must be in reality not only expansions in powers of a screening parameter but also 
the semiclassical ^-expansions, too. As stated earlier, we intend to restore the results 
of logarithmic perturbation theory for the nth eigenfunction and corresponding energy 
eigenvalue by means of the explicit expansions in powers of h. 

For this purpose we apply the method of series expansions in small parameter 
at the higher derivative, known from the theory of differential equations [^. Such an 
approach does not imply knowledge of the exact solution for initial approximation which 
is obtained automatically. Only the leading order in h for the energy eigenvalues need 
be determined first. With assumption that the energy eigenvalues and the logarithmic 
derivative of the wave function may be written as an asymptotic power series in the 
Planck constant, from (Q) we then have 

1 oo 

^ = A E Ekh'' (5) 

fc=0 
-| oo 

C{r) = -Y.C,ir)h"'. (6) 

Here we take into account that under substitution into equation(^) the coefficients of 
odd powers of h for the energy expansion (^) and even powers for (|^) are equal to zero. 

Before proceeding, a few words about the order in h of the combination he, 
appeared, as readily seen from (^, in our consideration. For velocities which are small 
in relation to the velocity of the light c, relativistic mechanics must go over into classical 
mechanics. So the non-relativistic energy E^r should be regarded as small quantity in 
comparison to the relativistic rest energy mc^. 
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Because for the Coulomb interaction, which is treated as the dominant part of the 
screened Coulomb potentials, the non-relativistic quantum-mechanical energy behaves 
as Enr ~ 1/^^, from the inequality Enr << mc^ we conclude that he ~ 0(1). It 
means that the factor he is merely the coefficient providing the right dimension. 

Therefore, to avoid additional complications, in foregoing consideration we put 
he = 1. It is worth noting that analogous statement for transition to classical limit is 
hold within the framework of the Dirac equation, too [|26| . 

Further, through the use of the ^-expansions (|) and @, on collecting coefficients 
of the same powers of h, from the Riccati equation (H) we obtain 



C'.ir) = - E', 

Co{r)Ci{r) = Eo[v{r) - E^] + mW{r) 

C[{r) + 2Co(r)C2(r) + Cl{r) = - E^ - 2E,E, + 2E,V{r) (7) 

+W'^{r) - V'^ir) 



Cfc-i(^) + E CAr)Ck-Ar) = - E E,E,., + 2E,.,V{r) k>2. 

j=0 j=0 

In the case of ground states, the recurrence system at hand coincides with one 
derived by means of the standard technique and can be solved straightforwardly. It 
will be recalled that complications of the logarithmic perturbation theory arise in the 
description of radial excitations when the nodes of wavefunctions are separated in some 
factor. We intend to circumvent this difficulty by making use of the quantization 
condition. Its fundamental idea that stems from the WKB-approach [|19], ^ is well 
known as the principle of argument in the analysis of complex variables. Being applied 
to the logarithmic derivative of wave function it means that 

j C{r)dr = hn ra = 0,1,2,... (8) 

where n is a number of nodes of the wave function, -R(r), and a contour of integration 
encloses only these nodes. 

This condition is exact and is widely used for deriving the high-order corrections to 



the WKB-approximation |]20| and the 1/ A^-expansions 0, ^. Now one important detail 
must be noted. The radial and orbital quantum numbers, n and Z, correspondingly, are 
specific quantum notions and need be defined before going over from quantum mechanics 
to classical physics. Therefore the quantization condition (H) must be supplemented 
with the rule of achieving a classical limit that stipulates the type of semiclassical 
approximation. 

It should be stressed that in our approach the quantization condition is quite 
distinguished from the WKB condition (|^) in two points: we apply the different rule of 
achieving a classical limit and choose another contour of integration in complex plain. 

Firstly, about the rule of achieving a classical limit. 
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Unfortunately, under the influence of standard textbooks on quantum mechanics, 
the semiclassical /i-expansions are usually associated solely with the WKB-approach, 
for which the rule of achieving a classical limit is 

h —>■ 0, n —>■ oo, / oo, hn = const, h I = const. (9) 

At that time the semiclassical 1/A^-expansion, being complementary to the WKB- 
method, requires the conditions ||^, |^ 

h ^ 0, n = const, / — > oo, hn 0, hi = const. (10) 

As it will be shown, the semiclassical treatment of logarithmic perturbative theory 
proved to involve the alternative possibility: 

^ — >■ 0, n = const, / = const, hn ^ 0, hi —>■ 0. (11) 

We should remark that the part of this rule concerned the orbital quantum number, 
/, differs from one used within the WKB-approach and the 1 / A^-expansion method. In 
our consideration this part, implying the first order in h for the quantity hi, has been 
used in deriving the system (|^). 

The second distinction is concerned with the choice of the path of integration in 
complex plane. In contrast to the WKB-approach and the 1/A^-method, we choose such 
a contour of integration which encloses not only the nodes of the wave function but the 
boundary point, r = 0, as well, and no other singularities. 

As it follows from equation (^, the regular branch of the function R{r) behaves in 
a neighborhood of the origin as 



R{r)r^o ~ r'/'+VK'-y.?+(^+m'u{r) , (12) 

where U{r) contains only nodes of the wave function. Hence, the direct integration 
of the logarithmic derivative R'{r)/R{r) around the origin yields the quantity 1/2 + 
\JWo — Vq + (/ + 1/2)2 addition to nodes in the quantization condition (|) which 
can be now rewritten as 

^ ^ C{r) dr = h(n + 1/2 + ^W^ - Vq^ + (/ + 1/2)^^ n = 0, 1, 2.. .(13) 

Taking into account that due to the rule ( pT]) the rihgt-hand side of this equality has 
the first order in h, and on substituting the expansion (H), our quantization condition 
takes its final form 

^jQ{r)dr= (n + 1/2 + ^W^ - V^' + (/ + 1/2)^) 6,^, n = 0, 1, 2...(14) 

where the Kronecker delta 6ij is used and the quantity, Wq — V^^ + (^ + 1/2)^, under 
square root sign should be not negative from the constraint to obtain a regular bound 
state wavefunction. 

A further application of the theorem of residues to the explicit form of functions 
Ci{r) easily solves the problem of the description of both ground and radially exited 
states. 
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3. Recursion formulae for perturbation expansions 

We proceed now to deriving the recursion relations for obtaining the nth eigenfunction 
and corresponding energy eigenvalue. 

Let us consider the system (0) and investigate the behaviour of the functions Cfc(r). 
From the first equation we have 



Co{r) = -^m^ - Ei , (15) 

where the minus sign is chosen from boundary condition. Then the function Ci(r) has 
a simple pole at the origin, owing to the Coulombic behaviour incoming in one or/and 
other part of the potential at this point, while the function Cfc(r) has a pole of the order 
k. Hence Cfe(r) can be represented by the Laurent series 

oo 

Ckir)=r-''Y.Cy , k>l. (16) 

With definition of residues, this expansion permits us to express the quantization 
conditions (0) explicitly in terms of the coefficients as 

Cl+' = N6o,k n = 0, 1, 2... (17) 



where N = n + 1/2 + ^Wi - Vq^ + (/ + 1/2)^. 

Thereby the common consideration of the ground and excited states has been indeed 
proved to be possible. 

The substitution of the series (|T6|) into the system (|^ and collecting coefficients of 
the like powers of r leads to the recursion relation in terms of the Laurent coefficients, 



2C^ 



p 

j=l p=0 



'0 
k 

+ Sk,i ^ EjEk-j — 2Ek-iVi-k+i 

j=0 

+ Sk,2 j2iKy^-p - WpW^-p) - S,,oSk,2 l{l + 1) (18) 

- 5k,i2mWi 

where for universality of designations we put Cq = Co(r), C° = 0, i > 0. 

In the case i ^ k, this formula is intended for obtaining C^, whereas if i = k, by 
equating the explicit expression for C^^^ to the quantization condition (|17D we arrive 
at the recursion formulae for the energy eigenvalues 



m 



r^Qr^l ( 1 rfc— 1 k 



E, 



k 



2{E,Cl + VoC^^)\Cl 



k—p 

.j=2 p=0 
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k 

+ 2Q{k - 2) ^ + 5k,MVr - W^W^) 

p=i 



(19) 



1 



j=l p=0 j=l 
2 

2E,_iri + 4,2 E(^.-^2-, - w^,w^2-,; 

j=0 



- bk \lmWx 



k>0 



where we use the step function 

e{k) =1, k>Q 
= 0, k<0 

and the plus-sign for Eq is chosen since the particle states enter to the bound state gap 
from the upper continuum. 

Thus, the problem of finding the energy-eigenvalues and eigenfunctions for the 
bound-state problem within the framework of the radial Klein-Gordon equation with 
screened Coulomb potentials can be considered solved. The equations ([TsD and (|l^) have 
the same simple form both for ground and exited states and define a useful procedure of 
successive calculation of higher perturbation theory orders convenient both for analytic 
evaluation and numerical computation by a computer. 

4. Discussion and examples 

In this section we consider specific application of our technique to perturbed hydrogenic 
problems. But first of all, in the proof of correctness of the recurrent formulae obtained, 
we check whether the exact results for the pure Coulomb potentials are restored. From 
equations (|19D it is readily seen that Eq does be exact solution to the Klein-Gordon 
equation with time- component of Lorentz four-vector and Lorentz-scalar Coulomb 
potentials |2^, [l| . As an additional check, we demonstrate that in this case our technique 
restores the exact solution for the eigenfunctions, too. 

Changing the variable from r to p, with p = 2r [m^ — E'^Y^'^ = 2/ir, and putting 
for simplicity, 7 = \J {I + 1/2)^ + Wq — Vq^ h = c = 1, from the system ([18|) we find 
that Cfc(p) = dkP~^, where 

dn = — 
" 2 

di=n+^ + j 

d2=nin + 27) (20) 



k-l 



dk = {1 - k) dk-i + E djdk-j k> 2 
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Integration of the functions Co(p) and Ci(p) gives exactly the part of the 
unnormahzed radial wave function 

Ri{p)=e-''/'p'/'+^Pn{p) (21) 

which provides the regular behaviour at the origin and correct one at infinity. 

The remaining part is a polynomial, P„(p) = J2k=o (^kP^, that satisfies the equation 

oo /I \ 

PLip)/Pn{p) = E d,p-' - (2 + 7) P-' (22) 

and whose coefficients are determined by the system 

00 

ak{n- k) + E CLjdj-k+i = . (23) 
j=k+i 

The combination of this equations, multiplied by a suitable dj with a view to taking 
into account equations (pOD, arrives at the following relation between two consecutive 
coefficients 

Qm^i ^ (m+ l)(m + 27) 
Um m — n — 1 

that is the recursion formula for the associated Laguerre polynomials L'^{p) (see, for 
example, |2^). Thus, described technique restores the exact result for the wave functions 
1^, too. 

Now, as an example of specific application of our technique we consider energy 
eigenvalues for the attractive Hulthen potential, having not only time-component of the 
Lorentz-vector, V{r) = — p^^y, but the Lorentz-scalar term, W{r) = —j^zj^ cis well. 
Because this screened Coulomb potential admits of the exact solution only in the case 



of s-wave p9| , pO|] , for describing orbital excitations we can use the perturbation theory. 
Then the coefficients of expansions (|I]) determining potentials are 

Vq = -a 



(25) 



Wo = -b 

and with applying equation (^) the analytic expressions for perturbation corrections 
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to the bound state energy take the form 

^ (Ny/m + a'^-h^-ah 



Eq 

El 
E2 

E4 



X 



2 m 



(am + bEo) 
{bm + aEo) 



24 (am + bEo) jj^^ni 
-E2 



?,{am + bEof - ^? L 



2 m 



{bm + aEo) 



X 



E. 



5760 {am + bEof fi^ 
{45 [{am + bEof + bb^ii'^{am + bEQf 

+L [24 a^m%Eo (6 E^ - m^) + 6 aW (tEq^ - 2 

-amn^{am + 2bEo) {5LE^ + 13Lm^ -Gm^) 

+62m (5 L m^ - 23 L + 6 E^) 

-6 6^ (a m + 6 Eof (l3 - 28 m^E^ + 5 

+6 b^E', (22 a^m^E^ - 2 a'm' - 5 fe^E^^)] } 

X^b (bm + aEo) 

~ 3840(am + 6Eo)V^^^ 

X {45 (am + 6^0)^ + 105 6V(a^ + &^o)^ 

+L [24 a^m^Eo (6 - m^) + 6 aW (tEq^ - 
—amjj^^am + 2b Eq) (5LEq + 13 Lm^ — 6m' 



2 m' 



+b^m^fi'^ (bLm"^ -23 LE^ + 6 El 
-2b\am + bEo)^ {l9 m'^ - U m'^E^ - S^o) 
+6 52^2 (22 a^m^E^ - 2 - 5 b'^E^)] } 



where /z^ 



m 



1/2)^ 



(26) 



£;2 , L = / (/ + 1), iV = n + 1/2 + y/H^o' - + 
It is readily seen that the use of ^-expansion technique leads to the explicit 
perturbation expansion in powers of the small parameter A. 

Besides, it is known the s-wave solution to the Klein-Gordon equation with vector 
and scalar Hulthen-type potentials, obtained in terms of hypergeometrical functions. 
From equation for allowed energy values, listed in [Q, one can derive the explicit 
expression for exact s-wave energy eigenvalues 



E 



-abm/n"^ + XaK^/2 + NKc^Jm'^ + Xmb - X^k^H'^/A 
N^ + ay{hc)^ 



(27) 



where = n + l/2 + ^byikc)'^ - a^/{hc)'^ + 1/4 and k = ^A^^ + a^l{nc)'^ - byikc)^. 

With performing the perturbation expansion of this expression in powers of A we 
arrive at the system (with he = 1) 

m —a b + N K 

En = — o = 

h^ + a2 
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^ _ X K{Ka + N b) 
^ ~ 2 iV2 + a2 

^ ^ Nk in' + &2) 

^ 8m iV2 + a2 ^ 
A /i^ 6 ^ 

-£^3 = 7^ E2 

2 m 

lb m^^ 

, b{3K' + 7b') 
^^--^ ^ 32 m3 



which coincides with ( pq) under the condition / = 0. Note that for s-waves there exists 
a critical value of the screened parameter A . From ( p7D the obvious condition 

m^c^ + Xmb- ^JIJL > (29) 
4 - 

yields 

2mc^ , , 

A,, = , . (30) 

\/N^h^c^ + a^ -b 
Then bound states can occur provided the condition A < Acr- 

In order to assess the speed and accuracy of the perturbation technique for the 
Hulthen potential with various ratios of its components we consider energy eigenvalues 
for the pure vector case, V{r) = —j^^; the pure scalar potential, W{r) = —j^^; 
and the equally mixed interaction. Typical results of calculation are presented in 
Tables where sums of first terms from our expansion for the energy eigenvalues for 
the bound state problem within the framework of the radial Klein-Gordon equation is 
compared with the results obtained by numerical integration, -Enum (in relativistic units 
h = c = m = 1). As the numerical integration procedure was used, the improved 



shooting method with the Noumerov integration scheme, described in literature |31 



which ensures the necessary exactness of calculation for the Sturm-Liouville problems. 

The Table 1 demonstrates dependence of accuracy of the perturbation expansions 
on the value of the screening parameter A. In Table 1 are shown the results of the 
calculations of the sum of first five terms from perturbation expansions for energy 
eigenvalues, computed for the first excited p-state {n = 1, I = 1) with parameters 
a = 1.0, b = 1.0 and the value of A varying from 0.05 to 0.15. It is seen that the 
accuracy of the perturbation description decreases with increasing the parameter A. 

The Table 2 illustrates the speed of the perturbation technique on the energy 
eigenvalue calculation for the Hulthen potential with various ratios of its component. The 
sequences of the partial sums of k corrections to the energy eigenvalues were computed 
for the states with I = 1, n = 1 and n = 2, and parameters a = 1.0, b = 1.0, A = 0.05. 

As can be seen from Table 2 the sequences of partial sums of perturbation 
corrections to the energy eigenvalues for the Hulthen potential have different behaviour 
in the pure vector case and for scalar case and mixed interaction. 
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Because the parameter b, which determines the scalar interaction, is involved 
linearly in all perturbation corrections of odd orders, starting from third (see equations 
(p6|)), these corrections in the pure vector case become equal to zero. Then we have one 
subsequence that tends smoothly to the exact value. 

When the pure scalar or mixed interaction is considered, we have two subsequences 
bounded below and above the energy eigenvalues. The average of these subsequences 
at the point of their maximal drawing together is proved to result in a quite good 
approximation to the exact value. 

5. Summary 

A new useful semiclassical technique for deriving results of the logarithmic perturbation 
theory for the bound-state problem withing the framework of the radial Klein-Gordon 
equation with the screened Coulomb potential having both time-component of Lorentz- 
vector and a Lorentz-scalar term has been developed. Based upon the ^-expansions and 
suitable quantization conditions, new handy recursion relations have been obtained. 
Avoiding the disadvantages of the standard approach these formulae have the same 
simple form both for ground and exited states and provide, in principle, the calculation 
of the perturbation corrections up to an arbitrary order in the analytic or numerical 
form. And, at last, this approach does not imply knowledge of the exact solution for 
zero approximation, which is derived automatically. As an example of application, 
perturbation expansions for the energy eigenvalues for the Hulten potential containing 
the vector part as well as the scalar component have been investigated with proposed 
technique. The different behaviour of the perturbation series in the case of potential in 
the form of time-component of a Lorentz- vector and for a Lorentz-scalar interaction has 
been found. 
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Table 1. The partial sums of 5 corrections to the energy eigenvalues Ey for pure 
vector Hulthen potential, Ew for pure scalar potential, and Ey+w for equally mixed 
interaction, and percentage error e, calculated with parameters a = 1.0, 6 = 1.0, in 
units U = m = c = 1. 



A 






Ey 




ev% 




Ew 






Ev+w 


ev+w% 


0. 


,05 


0, 


.95706870 


0, 


.00001 


0, 


.97392119 





.00003 


0, 


.84245453 


0.00001 


0. 


,06 


0, 


.96113964 


0, 


.00002 


0, 


.97739507 





.00008 


0, 


.85034939 


0.00002 


0. 


,07 


0, 


.96503895 


0, 


.00005 


0, 


.98064017 





.00020 


0, 


.85805052 


0.00004 


0. 


,08 


0, 


.96876527 





.00011 


0, 


.98365749 





.00045 


0, 


.86555969 


0.00008 


0. 


,09 


0, 


.97231703 


0, 


.00023 


0, 


.98644772 





.00091 


0, 


.87287852 


0.00016 


0. 


,1 


0, 


.97569245 


0, 


.00044 


0, 


.98901136 





.00171 


0, 


.88000851 


0.00030 


0. 


,11 


0, 


.97888956 


0, 


.00080 


0, 


.99134866 





.00307 


0, 


.88695104 


0.00052 


0. 


,12 


0, 


.98190615 


0, 


.00138 


0, 


.99345974 





.00527 


0, 


,89370738 


0.00086 


0. 


13 


0, 


.98473983 


0, 


.00229 


0, 


.99534457 


0, 


.00873 


0, 


,90027874 


0.00137 


0. 


,14 


0, 


.98738798 


0, 


.00369 


0, 


.99700303 





.01414 


0, 


.90666626 


0.00210 


0. 


,15 


0, 


.98984779 


0, 


.00577 


0, 


.99843492 





.02255 


0, 


.91287103 


0.00312 



Table 2. The sequences of the partial sums of k corrections to the energy eigenvalues 
Ev for pure vector Hulthen potential, Ew for pure scalar potential, and Ev+w for 
equally mixed interaction, calculated with parameters a = 1.0, b = 1.0, A = 0.05, in 
units h = m = c = 1. 



k n =1,1 = 1 n = 2,l = l 







Ev 




Ew 




Ev+w 




Ev 




Ew 




Ev+w 





0, 


.9341723590 


0. 


,9530618622 


0, 


.8000000000 


0, 


.9638612635 


0, 


.9726183555 


0, 


.8823529412 


1 


0. 


.9591723590 


0. 


,9768884088 


0, 


.8450000000 


0, 


.9888612635 


0, 


.9969338144 


0, 


.9294117647 


2 


0. 


.9570741392 


0. 


,9738581555 


0, 


.8423958333 


0, 


.9848180151 


0, 


.9915208448 


0, 


.9246200980 


3 


0, 


,9570741392 


0. 


,9739339118 


0, 


.8424609375 


0, 


,9848180151 


0, 


.9916561690 


0, 


,9247398897 


4 


0, 


,9570686998 


0. 


,9739202644 


0, 


,8424540955 


0, 


,9847983143 


0, 


.9916167626 


0, 


,9247202876 


5 


0, 


,9570686998 


0. 


,9739211933 


0, 


.8424545273 


0, 


.9847983143 


0, 


.9916195489 


0, 


.9247216080 


6 


0, 


.9570686381 


0. 


,9739209078 


0, 


.8424544790 


0, 


.9847977509 


0, 


.9916177074 


0, 


.9247213732 


7 


0, 


,9570686381 


0. 


,9739209397 


0, 


,8424544833 


0, 


,9847977509 


0, 


.9916179261 


0, 


,9247213972 


8 


0, 


,9570686368 


0. 


,9739209276 


0, 


.8424544828 


0, 


.9847977150 


0, 


.9916177253 


0, 


,9247213921 


9 


0, 


.9570686368 


0. 


,9739209295 


0, 


.8424544828 


0, 


.9847977150 


0, 


.9916177586 


0, 


.9247213928 


10 


0. 


.9570686367 


0. 


,9739209288 


0, 


.8424544828 


0, 


.9847977119 


0, 


.9916177282 


0, 


.9247213926 


Enum 


0. 


.9570686367 


0. 


,9739209289 


0, 


.8424544828 


0, 


.9847977115 


0, 


.9916177295 


0. 


.9247213926 



